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Abstract 



Exclusive hadronic S-meson decays into two-body final states consisting of a light 
pseudoscalar or vector meson along with an r/ or r]' meson are of great phenomeno- 
logical interest. Their theoretical analysis involves decay mechanisms that are 
^ ' unique to flavor-singlet states, such as their coupling to gluons or their "intrin- 

■ sic charm" content. These issues are studied systematically in the context of 

QCD factorization and the heavy-quark expansion. Theory can account for the 
experimental data on the B K^*^r]^'^ branching fractions, albeit within large 
uncertainties. 



1 Introduction 



The unexpectedly large branching fraction of order 6 ■ 10~^ for the inclusive production 
of high-momentum rj' mesons in B decays first reported by the CLEO Collaboration in 
1997 [H has triggered some theoretical activity aiming at explaining the dynamical origin 
of this enhancement. The large mass of the r]' meson as compared to other pseudoscalar 
mesons is one of the clearest manifestations of the axial anomaly in QCD. It is therefore 
natural to ask how the anomaly affects the production process of rj' mesons 0, ^, || . Other 
possible enhancement mechanisms (not mentioning those that invoke new fundamental 
interactions) include a large cc content of the 77' ^. Given the theoretical uncertainties 
involved in the calculation of semi-inclusive branching fractions, the dynamical details 
remain unclear except that there appears to be no fundamental problem to account for 
the observed branching fraction. 

Recent data on two-body final states containing rj' (summarized in Table |l|) show 
that the branching fraction for —>■ K~r]' decay is about six times larger than that for 
the corresponding B~ — > K~tt^ decay, confirming the semi-inclusive enhancement. The 
comparison with other two-body final states allows us to extract a more detailed pattern. 
For instance, if t]' is replaced by rj the corresponding branching fraction is suppressed 
rather than enhanced compared to the vr^ mode. Furthermore, the enhancement and 
suppression appears to be absent or even reversed when the pseudoscalar kaon is replaced 
by the corresponding vector meson, and with the current limited data no particularly 
striking pattern is visible when the kaons are replaced by pions or p mesons. 

The theoretical description of exclusive two-body decays has also improved, since the 
concept of factorization is now better understood [|1^]. The QCD factorization approach, 
in particular, has successfully explained the magnitude of tree and penguin amplitudes of 
final states with pions and kaons [O], while the calculation of strong interaction phases 
remains to be validated by experiments. It is therefore clearly interesting to see whether 
this approach can explain the pattern of branching fractions containing i]^'^ mesons as 
described above. The new element that has to be understood for this purpose is a 
fiavor-singlet amplitude, defined as the amplitude for producing a quark-antiquark pair 
not containing the spectator quark in the coherent fiavor state {uu + dd + ss) or a pair 
of gluons, where the quarks or gluons have small relative transverse momentum and 
hadronize into an r^^'^ meson. 

The analysis of this singlet amplitude, which has not been considered systematically 
in the QCD factorization approach so far, is the main purpose of this paper. We will 
indeed see below that the standard QCD factorization formula does not hold, but that a 
suitable modification allows us to obtain the fiavor-singlet amplitude at leading order in 
an expansion in l/m^ at the price of introducing one new non-perturbative parameter. 
We shall find that the QCD factorization approach can qualitatively account for the 
pattern of exclusive branching fractions described above, including the large rate for 
the Kt)' final state, but the theoretical uncertainties are rather large. In particular, 
we find that it is the constructive or destructive interference of non-singlet penguin 
amplitudes rather than an enhanced singlet penguin amplitude which is responsible for 
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Table 1: CP-averaged experimental branching ratios (in units of 10~^) on 
charmless B decays into two-body final states containing rj, rj' or 7r°. Upper 
limits are at 90% confidence level. 



Mode 



CLEO i 



BaBar 18 



Belle I 



B- 
B' 
B- 
B' 
B- 
B^ 



> K-r]' 
■ K-T] 



K 7r 

TT 







80 
89 



-10 , 7 

-18 
-16 



±9 



< 6.9 



< 9.3 



11.6 



-1-3.0-1-1.4 
-2.7-1.3 
14 fi+5.9+2.4 
-■-^■"-5.1-3.3 



67±5±5 
46 ± 6 ± 4 



12.8 ± 1.2 ±1.0 
10.4 ± 1.5 ±0.8 



77 Q 4-6.2 +9.3 
' ' -^-5.9-8.7 



68.0 



-10.4+8.8 
- 9.6-8.2 



5 2+i-'^ 

^"^-1.5 



13.0+1^ ± 1.3 
8.0^|-?±1.6 



B- 
B^ 

B- 
B' 
B- 
B' 



> K*-r]' 
' K*-r] 
■ K*-n^ 



< 35 

< 24 
26.4+^:^ ±3.3 
13.8^^1 ± 1.6 

< 31 

< 3.6 



< 13 
< 33.9 
19.8+^1 ± 1.7 



26.5 



-7.8 
-7.0 



±3.0 



16.5+^2 ± 1-2 



the distinctive pattern of strangeness-changing transitions to r/*^') mesons, much along the 

The improved description of data compared 
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lines envisaged qualitatively by Lipkin 
to the naive factorization analysis [|T^ comes from the radiative enhancement of the 
non-singlet penguin amplitude, which also underlies the sizeable branching fractions 
for TiK final states. Our findings are in contrast to other recent analyses of the rf^'^ 
modes |T^ , |15|, |T6|, where the large branching fractions for the Kr]'^'^ final states were 
attributed due to an enhanced fiavor-singlet penguin amplitude. Similarly, the singlet 
amplitudes we obtain are smaller than those inferred from phenomenological analyses 
using SU(3) symmetry and certain dynamical assumptions about fiavor topologies [|^. 
For completeness we note that the decays B Kt]^'^ have also been analyzed using the 
perturbative QCD approach [1^, however no singlet-specific mechanisms besides rj-rj' 



mixing were investigated in this study. 



2 Implementation of rj—r]' mixing 

In the calculation of weak decay amplitudes with an i]^'^ meson in the final state, we need 
several matrix elements of local operators evaluated between the vacuum and r]^'^ . These 
are the matrix elements of the flavor- diagonal axial-vector and pseudoscalar current 
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densities, 



{P{q)\qYim\^) = -A^fpq^, 2m,{P{q)\qj,q\0) = -^h%, 

{P{q)\-si^l,s\Q) = -if'pq^, 2m,(P(g)|s75s|0) = -ih'p , 

where q = u or d. We assume exact isospin symmetry and identify = ^{rriu + m^). 
We also need the anomaly matrix elements 

{P{q)\^G^^G^nO) = ap, (2) 

where we use the convention 

1 
2 



for the dual field-strength tensor. In all cases P = t] or t]' denotes the physical pseu- 
doscalar meson state. We also need the generalization of the local quark operators to 
the corresponding light-cone operators, which define the twist-2 and twist-3 light-cone 
distribution amplitudes. The treatment of mixing for these generalizations will be anal- 
ogous to the case of the local operators. Finally, we need the current matrix elements 
{P\qrb\B) , which we decompose into Lorentz-invariant form factors as for any other 
pseudoscalar meson. 

The equations above define ten non-perturbative parameters fp, h'^p, and ap, which 
however are not all independent. Taking the divergence of the flavor-diagonal axial- vector 
current, 

dMl'im) = 2tm, q^,q - ^ G^, G^'>^^ (4) 

(and similarly with q replaced by s) yields four relations between the various parameters, 
which can be summarized as 

ap= ^^~J^^ =/^l.-/^m^. (5) 

Without further assumptions this leaves us with six independent parameters. It is con- 
ventional to write each set of two parameters corresponding to P = r], t]', such as {/^, /^,}, 
in terms of a parameter (such as fs) and a mixing angle (such as 0^). A general treatment 
then implies three such parameters and three independent mixing angles. 

In the SU(3) flavor-symmetry limit, where \7]) is a flavor-octet and \ri') a flavor- 
singlet, it follows that = -V2f^ and f^, = f^,/V2, /i^ = -v^/i? and h'^, = /i^,/\/2, 
and = 0. However, it is known empirically that SU(3)-breaking corrections to these 
relations are large. In the following we shall not rely on SU(3) flavor symmetry but 
instead introduce another assumption, expected to be accurate at the 10% level, to 
reduce the number of hadronic parameters. This assumption leads to what will be 
referred to as the Feldmann-Kroll-Stech (FKS) mixing scheme |]19[. In the absence of 
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the axial U(l) anomaly, the flavor states \r]q) = {\uu) + \dd))/\/2 and \r]s) = \ss) mix 
only through OZI-violating effects known phenomenologically to be small. We therefore 
assume that the anomaly is the only effect that mixes the two flavor states. (This 
assumption implies, in particular, that the vector mesons uj and are pure {uu + dd) 
and ss states, respectively, as is indeed the case to very good approximation.) In a chiral 
Lagrangian treatment of the pseudoscalar mesons (including the r]' meson, which can be 
done in a combined chiral and l/N^ expansion [^), the anomaly introduces an effective 
mass term for the system oir]'^'^ states that is not diagonal in the flavor basis {\r]q), |r/s)}, 
and since this is by assumption the only mixing effect, the FKS scheme amounts to a 
scheme with a single mixing angle in the flavor basis. If the physical states are related 
to the flavor states by 

(6) 

then the same mixing angle applies to the decay constants fp and h^p with the normal- 
ization given by (|l]). We therefore write 

f^ = fqCos(f), /' = -/,sin0, 

(7) 

and an analogous set of equations for the hp. This defines four new parameters fg^s 
and hq^g. Inserting these results into (|^) allows us to express all ten non-perturbative 
parameters in terms of the decay constants fq, fg and the mixing angle (f). We obtain 

hq = fq [rn^n cos^ + m^, sin^ 0) — v^/s (m^, — m^) sin cos , 

f (8) 
72 



hs = fs {fTT'l' cos^ + sin^ 0) ^ (m^, — m^) sin cos ( 



and 



ttri = ifg ''^ri ~ ^q) COS = ^= — - sin cos {—fq sin + V2fs cos ( 

1 , „ o . . 'ml' - fnl 



(9) 



a^/ = -= {fq m'^, — hq) sin = ^^^~~/= — ~ ^ "^^^ *^ (fq '^'^^ ^ + V^/s sin 0) . 

V 2 V 2 

The three remaining parameters of the FKS scheme have been determined from a fit 



to experimental data, yielding |19 



/, = (1.07 ±0.02)/,, /, = (1.34 ±0.06) A, = 39.3° ± 1.0°, (10) 

where the errors do not include a possible systematic uncertainty from the theoretical 
assumptions underlying the FKS scheme. The numerical values of the other parameters 



4 



Table 2: Decay constants (in MeV), pseudoscalar densities and anomaly matrix 
elements (in GeV^) of r^^'^ in the FKS mixing scheme. 



hq 


0.0015 ±0.004 


fi 


108 ±3 


K 


0.001 ±0.003 




0.087 ±0.006 


Jr] 


-111 ±6 


K 


-0.055 ±0.003 




-0.022 ±0.002 


Jri' 


89 ±3 




0.001 ±0.002 




-0.057 ±0.002 


frt' 


136 ±6 


11^, 


0.068 ±0.005 



derived from these inputs are given in Table 0. As can be seen from the table the 
expression for hq in (|^) is not useful to compute this parameter accurately in practice, 
since large cancellations occur between the terms proportional to fq and fs- In fact, from 
(HI) it follows that hq must vanish in the limit and so is expected to be very small 

{hq/hg ^ mq/ms).^ An alternative expression is obtained with the help of the chiral and 



large-Nc expansion. To leading order this gives [20 



fq = U, /. = ^2/i - /2 = 1.41A , ^^^^ 

hq = fqml = 0.0025 GeV^ , h, = f, {2m\ - ml) = 0.086 GeV^ . 

These results are also prone to large errors, since the chiral expansion is expected to 
apply only if the masses of the pseudoscalar mesons (now including the r]') are much 
smaller than the masses of the other hadrons. We conclude that hq is poorly determined 
at present. Note that although hq is small it cannot be neglected in general, since some 
decay amplitudes are proportional to hq/rriq. 



3 The flavor- singlet amplitude 

A distinctive feature of hadronic B decays with an rj^'^ meson in the final state is that 
these states contain a flavor-singlet component in their wave function, which opens pos- 
sibilities for novel decay mechanisms that do not occur if the final state contains only 
flavor non-singlet states such as pions and kaons. In this section we discuss the different 
contributions to the flavor-singlet amplitude in the QCD factorization approach. 

Neglecting weak annihilation terms for now, we write the B — >■ KP decay amplitudes 
(where P = ri^''>) as 

A{B ^KP) = ^^Y1 MKP) , (12) 

* p=u,c 

"'^Because the masses of and r]' are themselves functions of the mixing angle it would be inconsis- 
tent to interpret (||) as equations for hq and hs as functions of (j) with and m,,' fixed to their physical 
values. Rather, only those values of m^, m^/ and (j) leading to a small value of hq are acceptable. Setting 
hq ={) gives (j) ~ 38°, which is close to the phenomenological value quoted above. 
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where Xp^ = VpbV*g, and 



+ fp [5pc a2{KP,) + aliKP,)] | . (13) 

Here Fo^^*^(0) is a P — >■ M transition form factor evaluated at = 0, qg = u or d 
denotes the flavor of the spectator antiquark in the B meson, Cu = 2/3 and = —1/3 
are the corresponding electric charge factors, and the symbol 6q^u (and similarly 6pu) is 
defined such that 6q^u = 1 ii qs = u and zero otherwise. The notation Pi will be explained 
below. 

The various coefficients in the expression for the amplitudes correspond to different 
flavor topologies contributing to a given decay The arguments in parentheses are 
such that the first meson inherits the spectator antiquark from the B meson, while 
the second meson is the "emission particle" produced at the weak vertex. In naive 
factorization, the are given in terms of combinations of Wilson coefficients of the 
effective weak Hamiltonian as 



a,{PK) = C^ + j^^, al{PK) = C, + j^^+r'^ [C, + j^^ 



a2{KP,) = C2 + ^ , al{KPs) = C4 + ^ + <M ^6 + 

and 



(14) 



al{KP,) = C3 + ^ - + ^ ) , (15) 



where = 2m\/[mb{ms + rnq)\ in the definition of a^'^{PK), and r^" = / {ff mbrris) in 
the case of a^{KPs). Analogous expressions hold for the electroweak penguin coefficients 
'^s.EW Q^4EW Note that the singlet amplitude a^{KPi) is generated by penguin 
operators in the effective weak Hamiltonian through transitions of the type b — >■ sqq, 
where the gg-pair hadronizes into an r]^'^ meson. 

The calculation of corrections of order as to the coefficients ai,2,4 is analogous to the 
calculation for pions or kaons discussed in |[ll|]. (In the notation of that paper, we have 
ai,2 = CLi,2, «3 = 03 — as, 04 = a4 + r^ag, as.EW = og — clt, and q;4_ew = ^lo + fx^»-) 
In the FKS scheme, the rj^'^ system is described by two leading-twist quark-antiquark 
light-cone distribution amplitudes, (pq{x) and 0s (x), defined in analogy with the decay 
constants fg and fs- The subscript on Pg and Ps in (|T3|) means that the corresponding 
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Oi 

Figure 1: Two-gluon emission from a quark loop. A second diagram with the 
two gluons crossed is imphed. 



distribution amphtude must be used in the expressions for the ai parameters. The 
one-loop vertex corrections to the quark-antiquark amplitude and a contribution from 



spectator scattering to the singlet amplitude have also already been computed in 
This result is not complete, however, since it ignored the possibility that the rj'^'^ meson 
is formed from two gluons. 

In the following subsections we discuss three contributions to the singlet decay am- 
plitude related to the gluon content of the r]^'^: the h sgg amplitude and its relation to 
an effective charm decay constant, spectator scattering involving two gluons, and singlet 
weak annihilation, the latter being suppressed by at least one power of A/nib in the 
heavy-quark limit. Before we begin this discussion we comment on the treatment of the 
heavy-quark content of the t]^'^ meson. In QCD factorization the decay amplitude is 
factorized into short- distance kernels and light-cone distribution amplitudes. When the 
factorization scale is smaller than the heavy-quark mass, heavy-quark contributions are 
part of the short- distance kernels, and the light meson is described in terms of distribu- 
tion amplitudes of light quarks and gluons. If, in the limit rrif, oo, the charm-quark 
mass is held fixed, one should still introduce charm light-cone distribution amplitudes 
to sum large logarithms of mh/nic, but these distribution functions (and the associated 
decay constants) can be computed in terms of those of light quarks. In ([T3|) we included 
the possibility of such an induced charm contribution in the form of terms ai{KPc), 
to which we will return below. 



3.1 The b — > sgg amplitude 

The b sgg amplitude for general gluon momenta has been calculated in in the 
electroweak theory. What we need here is the special case where the two gluons have 
small invariant mass to form an r]^'\ and where the 6 — s transition is induced by an 
operator in the effective weak Hamiltonian. The relevant diagrams are shown in Figures |I| 
and ^ 

The diagram of Figure |l| is only relevant when the internal quark is a charm or bottom 
quark. (Top quarks have already been removed from the effective weak Hamiltonian.) 
Diagrams with light-quark loops are governed by hadronic scales and their short-distance 
part consists of the Oi insertion only. This is already taken into account in the naive 
factorization results (|1^ and (|T5|). For an operator O = {sTib){qT2q), where the Fj 
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denote arbitrary spinor and color matrices, we find (assuming the two gluons to be in a 
color-singlet configuration) 



A{b ^ sgg) iFig. I = -{usT.Ub) {g{qi)g{q2) |tr {T^ A) |0) 



(16) 



where Ui denote quark spinors, and the trace is performed over spin and color indices. 
The quantity A, which is proportional to the 3x3 identity matrix in color space, is given 
by 



A 



AixNr Sm„ 



[1 - 4r) F(r) + 4r 

1 r 



nr)^l,G%G^^>^'' 



+ 



1 - Fir) 



(17) 



where nig is the quark mass in the loop, r = rrig/q'^ — ie, and = (gi + ^2)^ will later be 
identified with the square of the pseudoscalar meson mass. The function 



Fir) 



Ar arctan 



V4r - 1 



1 



1 

12^ 



0{l/r' 



(18) 



describes the dependence on the invariant mass of the two gluons and the quark mass 
in the loop. The last term in ([T7| ) may be simplified using the identity 

(19) 



qA qA,ij, ^ 9a/3 qA qA,^j,u 



fia 



which for local operators is implied by the antisymmetry of the field-strength tensor. In 
the limit where 3> g^ (i.e., for g = c, 6), we obtain 



A 



AirNr 12m 



1 



- 6mg) 
48m2 



G^^,G^'^^ + 0{l/m^ 



(20) 



For a pseudoscalar state only the matrix element of G^^^ G^'^" is non-zero. Evaluating 
the Dirac and color structures for the various operators in the effective weak Hamiltonian, 
we obtain for the contribution to the B — > KP decay amplitudes 



A. 



charm 



ap 
I2ml 



C4 — Cfi 



(21) 

where q = p — p' and where we have neglected the very small contribution from bottom- 
quark loops. The definition (^ has been used to evaluate the matrix elements of the 
anomaly. It follows from (|^) that, to a very good approximation (i.e., neglecting the 
small quantity hg), ap/m?p = —fp/\^. Comparing (|2T|) with (p!^)-([T5|), we see that the 
charm-loop contribution can be incorporated into (|13D if we identify 



ff 



m 



2 fQ 

p Jp 



ap 

8 



(22) 



Oi O, Os, Os, 

Figure 2: Contributions to the h — > sgg amplitude that are power suppressed 
in the heavy-quark hmit. 

with the leading-order result (in an expansion in 1/ mc) for a new decay constant defined 
as 

(P(g)|c7^5c|0) = -^/^g^ (23) 

which might be interpreted as a measure of the "intrinsic charm" component of the rf^'^ 
wave function. With our definitions fp is negative and takes values ~ — IMeV and 
f^i —3 MeV. These results for the decay constants agree with a different derivation in 
A similar discussion of the charm-loop diagram in the context of non-leptonic B 



decays has already been given in 13 



Note that the charm decay constant is formally power-suppressed for heavy charm 
quarks (but not power suppressed if only mf, goes to infinity), but that the factor of has 
disappeared in the evaluation of the anomaly matrix element. The charm contributions 
are potentially non-negligible, since they multiply the large Wilson coefficients Ci^2- The 
eventual smallness of this contribution is due to the factor 1/12 in (|2^) , and because the 
large Wilson coefficients enter in the color-suppressed combination C2 + Ci/Nc- Since 
we have identified a new leading-order contribution to the factorized amplitude ([T3|), 
we should in principle consider corrections to a2^3{KPc), which would require the 
calculation of two-loop diagrams. We expect that these diagrams can be factored into 
genuine corrections and a contribution with the structure of a one-loop kernel (with 
external charm quarks) folded with a calculable charm-quark distribution amplitude. 
Due to the smallness of the leading-order charm contributions, however, the calculation 
of this correction is not required within the accuracy of our analysis. 

There are other diagrams that could potentially contribute to the b sgg amplitude, 
in which one or both gluons are emitted from external quark lines. Figure ^ shows the 
relevant graphs with one gluon emitted from a penguin loop or from the chromo-magnetic 
penguin operator, and the other gluon emitted from the b or s-quark lines. The loop 
diagrams are non-zero only for an insertion of the penguin operator O5. In order to 
obtain a leading-power contribution to the decay amplitude we must pick out the large 
components of all momenta in the above expression. However, setting qi = xq + . . . , 
q2 = xq + . . . and pb — Ps = q + ■ ■ ■ one finds that the leading term is proportional to 
thus vanishes for on-shell gluons. It follows that the diagrams in Figure ^ 
are power suppressed with respect to the diagram in Figure |I|. A similar analysis shows 
that also the graphs where both gluons are emitted from the b and s-quark lines yield a 
power-suppressed contribution. 



9 





Figure 3: Spectator-scattering contributions to the B K^*^?]^'^ decay ampli- 
tudes. The shaded blob represents the r]^'^g*g* form factor. 

3.2 Spectator mechanism 

The mechanisms described so far in this section provide the leading contributions to the 
flavor- singlet amplitude in the heavy-quark limit, which start at zeroth order in as (naive 
factorization). In practice, these contributions are numerically rather small, because 
they are suppressed by the small Wilson coefficients C3...6 or the mass ratio mp/(12m^). 
We will now identify another leading-power contribution to the fiavor-singlet amplitude, 
which starts at order but is free of such suppression factors. This contribution is 
related to the spectator-scattering diagrams shown in Figure ^ 



The fiavor-singlet spectator mechanism was considered first in |[14| , |I5|, |I6| using the 
perturbative QCD approach. In fact, the kinematics of the spectator graph is such that 
at leading power the two gluons cannot be coUinear with the outgoing i]^'^ meson, and 
at least one of them must be a (semi-) hard gluon. Specifically, assigning momentum 
Pq = yPK + k± + . . . to the antiquark in the kaon (with y = 1 —y) and a soft momentum 
/ to the spectator quark in the B meson, it follows that for generic values of y not 
close to or 1 the virtualities of the two gluons are ~ —'^VPk ■ ^ = 0(mf,A) and 
qf ^ ym% = 0{ml). The coupling of two hard gluons to the r]^'^ meson can be described 
by a perturbative form factor F{ql,q2) = 0{as) defined in terms of the r]^'^g*g* vertex 

as n,!! 

^Fp{ql, ql) e^,„^ ql^ 5ab (?i) 4*(?2) , (24) 

where P = rj' or rj. It follows that for generic values of y the hard spectator-scattering 
contribution is an effect of order and thus is beyond the accuracy of a next-to-leading 
order analysis. 

Nevertheless, it is instructive to work out the result for this contribution. In the 
kinematic region where \qf\ ^ Ig^l the leading contribution to the form factor can be 



written as 



where 



Fp{qlql) = -^,Cp, {\ql\:>\ql\) (25) 

Cp = V2fl fdx ^ + rP fdx ^ + . . . (26) 

Jq bxx Jq bxx 
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contains convolution integrals involving the leading-twist quark-antiquark distribution 
amplitudes of the t]'^'^ meson mentioned earlier (twist-3 quark-antiquark distribution 
amplitudes do not contribute due to chirality conservation), and the dots represent a 
contribution involving the two-gluon distribution amplitude, which we omit for simplic- 
ity. In the hmit of asymptotic distribution amplitudes this latter contribution vanishes, 
and Cp \/2fp + fp. At leading order in 1/mf,, we find for the hard spectator-scattering 
contribution to the B Ki]^'^ decay amplitude 



^spec 



Am 



Cp fpfK 



mp 



Uy 



B Jo y 



my)My) + r^Piiy)<l^?iy)] 



(27) 



where 4>K{y) ~ Qyy is the leading-twist distribution amplitude of the kaon, 4>p{y) ~ 1 
is one of the twist-3 amplitudes, and the hadronic parameter is defined as the first 
inverse moment of one of the two leading i?-meson distribution amplitudes |TD|, [TT[]. The 

OT*Ci -flno v^Qvi m Tim Ij'OT'm olo Mn^TT\rlTi r^tz^ri im (/l(Jl o tt t~\ (R/ll r\T IH 111 

we have 



quantities -Rf3(y) are the penguin kernels introduced in (49) and (54) of [|Tl|]. Explicitly, 



P!{y) = c, 

+ (^4 + Q 



4: rrib 2 ' 
o In — + - - G{sp,y} 

6 jJ, 6 



|ln^ + ^-G(0,y)-G(l,i/) 
3 /i 3 



1^ In - (n^ - 2) G(0, y) - G(.e, V) - G(l, y) 



2Q 



■off 



y 



(28) 



with Sq = {mq/mby, Uf = 5, and the penguin function G{s,x) as defined in the same 
reference. The expression for P^{y) is identical to that for -Pf (?/)) except that the term 
—2Cgg/y must be replaced with —2Cgg. Inspection of the above expressions shows that 
the perturbative hard-scattering contribution has a logarithmic endpoint singularity as 
?/ — > 1, corresponding to the region where the gluon exchanged with the spectator quark 
becomes a soft gluon. For the contribution proportional to CgJ both the leading and 
subleading- twist terms are divergent, whereas for the penguin loop contributions only 
the twist-3 term diverges. In []T3|, |TB[ these singularities were regularized by keeping 
subleading terms in the quark propagators, which effectively corresponds to adopting 
a model for the non-perturbative endpoint region. While this may be legitimate for 
obtaining a rough numerical estimate of the spectator-scattering effect, we do not sub- 
scribe to the conclusion reached by these authors, that the diagrams in Figure |^ are 
free of endpoint singularities and thus short-distance dominated. On the contrary, the 
logarithmic sensitivity to the endpoint region indicates that the soft spectator-scattering 
mechanism has one less power of associated with it, and so it is formally leading with 
respect to the hard spectator-scattering contribution. 

The appearance of a leading-power soft spectator-scattering contribution is a novel 
feature of our analysis, which is specific to the case of a light flavor- singlet final-state 
meson. The QCD factorization formula for B decays into two non-singlet mesons must 
be extended to account for this effect. We will now argue that factorization still holds in 
a generalized sense: the soft spectator-scattering contribution can be parameterized in 
terms of a convolution of the light-cone distribution amplitudes of the singlet meson with 
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a perturbative kernel (the expression Cp in (|26|)), multiplied with a B ^ K "form factor" 
defined as the matrix element of a non-local operator. For the purpose of illustration we 
focus on the leading-power twist-2 contribution, corresponding to the term proportional 
to Cg^ in (^). Momentum conservation implies that in the endpoint region (where the 
momentum q2 in Figure ^ is soft) the gluon emitted from the weak vertex is still semi- 
hard, since qf = (pp— q'2)^ ~ —'^'Pp'12 = O(mbA). In the limit where \q\\ ^ l^al = O(A^), 
the leading contribution to the i]^'^g*g* form factor is still given by (p5D, however now one 
of the two factors Qs is a non-perturbative coupling. Effectively, the soft gluon couples to 
a compact (qq) pair in a color-octet state, whose internal features it cannot resolve. The 
dependence of the soft spectator-scattering mechanism on the soft gluon momentum q2 
can be obtained from a weighted integral over a dual gluon field-strength tensor along 
the light-like trajectory of the r]^'\ using the formalism introduced in Section 4.3.2 of the 



second paper in |]I0|. By dimensional analysis, the resulting non-local matrix element 
scales like a heavy-to-light form factor at large recoil, and we thus parameterize it as 
(we suppress the Wilson lines required to render the non-local matrix element gauge 
invariant) 

2 

dss{K\rsiO) [^,7^](1 + ^,)gsn''G,o.{sn) b{0)\B) = !!^^/-^(0) , (29) 

where is a null-vector in the direction of the r]^'^ meson, and the "form factor" 
JF^^^(O) scales like {K/nihY^'^ the heavy-quark limit. With these definitions, the 
result for the soft spectator-scattering diagram can be represented as a contribution to 
the coefficients af(Pg,s) in ([131) , 

«3(^n,.)Utspec--^^Cg,(/i.)^ dx^— + ...j ^p^^(Q), (30) 



where = \/^bAh with Kh = 0.5 GeV serves as a typical scale for the semi-hard 
gluon propagator, and as above we have neglected a contribution from the leading-twist 
two-gluon distribution amplitude of the meson P. This result shows explicitly how 
the factorization formula must be modified to account for the effect of soft spectator 
scattering. 

A rough estimate for the form factor JF^~*^(0) can be obtained by regularizing the 
leading-power logarithmic endpoint divergence in (^7[) using a cutoff such that y > 
^h/^B- This yields 



(0)«^^^l"^-lj«0.3, (31) 

which is of about the same magnitude as the conventional form factor F^^^{0). The 
resulting contribution to a^{Pg^s) is then about 2 • 10^^, roughly half as big as the naive 
factorization contribution in (|15D . 
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Figure 4: Representative singlet annihilation diagrams. 



We remark here that we are unable to reproduce the large enhancement of the 
spectator-scattering effect for K mesons relative to K* , which according to [1^ is respon- 
sible for the large Krj' branching fractions. In this work the difference between pseu- 
doscalar and vector mesons arises when twist-3 distribution amplitudes are included. 
We find that including the term proportional to C|J in P^{y 
coefficient of the logarithm in (^) changes from 1 to 1 -|- 
mains unaltered for K*. We suspect that the huge effect observed in |16] is related to an 
error in the kernel for the twist-3 contribution, which leads to an artificial 1/y ~ tjib/ ^ 
enhancement of this contribution. 



into our estimate, the 
: 1.15 for kaons and re- 



3.3 Singlet weak annihilation 

The discussion of weak annihilation effects in flavor non-singlet non-leptonic decays has 
received considerable attention recently |ll], |2^, Annihilation effects belong to the 
class of power-suppressed corrections to the decay amplitude, which cannot be calculated 
in the factorization approach because they receive leading contributions from soft par- 
tons. The analysis of the flavor non-singlet penguin amplitude a\{MiM2) has shown that 
annihilation effects are probably not large, yet they constitute the largest uncertainty in 
the calculation of this amplitude |11|. We now consider the corresponding effect in the 



flavor- singlet case. 

Formally, the leading singlet annihilation amplitude is generated by the emission 
of two gluons, which form the rj^''' meson. Figure ^ shows three representatives of the 
twenty relevant Feynman diagrams. The six diagrams where both gluons couple to the 
constituents of the B meson are part of the two-gluon contribution to the B i]^'^ form 
factor used as an input in the factorization formula, and need not be considered further. 
Heavy-quark power counting of these diagrams shows that the leading term comes from 
the diagrams with two gluons emitted from the spectator quark, and confirms that this 
contribution is of the same order in the l/ruf, expansion as the quark contribution to 
the form factor. To compute the remaining terms in the hard-scattering approach, we 
define the leading-twist two-gluon light-cone distribution amplitude by (omitting again 
a path-ordered exponential) 



{P{q)\G/'^{w) Gj,(^)|0) = {V2fl + /^) q,q, dx e^^^^'^^"^-^ 0f (x) , (32) 
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where {z — wY = 0, and the particular combination \/2fp + fp is proportional to the 
flavor-singlet decay constant. The corresponding momentum-space projection onto the 
on-shell two-gluon scattering amplitude (with the gluon polarization vectors e* removed) 
reads 

lj;i.M,;.W,^,t^m, (33) 
2 — 1 q ■ n XX 

where q is the momentum of the r]^'^ meson, is a light-like vector in the direction of the 
kaon, and the index a {(3) belongs to the gluon with momentum xq (xq). In the FKS 
scheme we set 0^ (x) = (pg{x) = (f)g{x). The first term in the Gegenbauer expansion of 



x) reads |28, Eg, RO, m 



bg{x) = 5B^{fi)x'^x'^{x-x) + ... . (34) 



The coefficient -BK/^) as well as all higher Gegenbauer moments vanish for fi —>■ oo. The 
process 77* —> rj^'^ can in principle constrain i^K/^)- In our normalization convention^ 
the analysis performed in gives -Bf (1 GeV) = 2±3. Note that the second Gegenbauer 
moment of the singlet quark-antiquark amplitude is much smaller, of order —0.1, which 
underlines the weak sensitivity of the process to the gluon contribution. 

The result of the computation shows that we obtain a A/ nib correction to the singlet 
decay amplitude from the diagrams where one gluon couples to the spectator quark and 
the other to the constituents of the kaon. The remaining diagrams are further suppressed 
in 1/mfe. However, for the graphs where both gluons are emitted from the constituents 
of the kaon,0 this additional suppression is compensated by the "chiral enhancement 
factor" . For this configuration we also find that the result suffers from endpoint 
singularities. As in the flavor non-singlet case we can then only give an estimate of 
annihilation effects, based on cutting off the convolution integrals at parton momenta of 
order A. To simplify this estimate, we neglect electroweak penguin annihilation and note 
that annihilation through the tree operators Oi^2 competes with large tree amplitudes 
and is CKM-suppressed for 6 — > s transitions. We therefore concentrate on annihilation 
effects that contribute to the singlet penguin amplitude. In this approximation, these 
are accounted for by substituting 



where 



aliKPg,,) ^ aliKPg,,) + Pss , (35) 



. ,^,C,\ fBfxr^ 27ia, [\ ^(y) x - x 



Nj m%Fo^^^{0) io yy Jo " ' x^x^ 



^Thc definition implies that our distribution amplitude is ^ y/Cp/uf times the distribution 
amplitude defined in Appendix A of and Cp/i'^nf) — 2/9 times the distribution amplitude assumed 
in Sections 3 and 4 of the same paper, to which the determination of -B|(/i) refers. 

■^The singlet annihilation contribution where both gluons are emitted from the kaon can be interpreted 
in terms of a time-like Krj'^'^ form factor, evaluated at = rri^g. 
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and a small contribution proportional to C^/Nc has been neglected. Note that (5sz 
is formally power-suppressed by a factor r^/rrib ~ 1/ml relative to a^, similar to the 
suppression of the (numerically) dominant penguin annihilation amplitude for the flavor 
non-singlet case. Within the approximations described above, our result for (3s3 refers 
to penguin annihilation into a pseudoscalar meson {B K) with both gluons radi- 
ated from the constituents of the kaon, as indicated in the last diagram in Figure ^ 
The corresponding diagrams for annihilation into a vector meson {B K*) vanish for 
asymptotic distribution amplitudes. 

To obtain a numerical estimate of /J^s we insert ( p4D for the two-gluon distribution 
function, the asymptotic twist-3 distribution amplitude (f)p{y) = 1 for the kaon, and 
parameterize the endpoint- divergent integrals by = j} dy/y as in This gives 



^^JV~ Jo '^'^^^^^i^ " T ^ ^ 

Evaluating all quantities at the scale fih = y/rribAh with = 0.5 GeV, and choosing 
Xa = ln{mB/ Afi), we find 

Ps3 - -9 ■ 10-^ Bl . (38) 

The value of jSss is very uncertain. Using i?| = 2 ± 3, and allowing Xa to vary from to 
twice the estimate given above, we obtain —9 • 10~^ < j3s3 < 2 ■ 10^^, but the actual size 
may be much smaller, since the central value for the gluon Gegenbauer moment may not 
reflect the true magnitude of -Bf . 



4 Phenomenological implications 

We now discuss the phenomenological consequences of our results, present numerical 
results for the various branching fractions, and investigate to which extent the striking 
patterns seen in the experimental data in Table ^ can be understood in the context of 
QCD factorization. 

The leading contributions for the B — >■ Krj^'^ decay amplitudes have already been 
given in (|T2p and (|TB|). The corresponding result for the case where P = reads 



Ik / 

71V 



2 



+ m|Fo^-^(0) ^ !^6,^a2{K7r) + ^al^^{K7r)^ , (39) 

where Qs = u or d denotes the flavor of the 5-meson spectator quark, and = 1 
and (Tc; = — 1 are sign factors. At subleading order in the heavy-quark expansion the 

"^By dropping the contribution proportional to C3 we eliminate the diagrams with a gluon coupling 
to the _B-meson spectator quark, which we cannot estimate in the conventional way, since the hard scat- 
tering amplitude depends on both the plus and minus components of the spectator-quark momentum, 
rendering the light-cone projection onto the B meson as defined in p3|, M invalid. 
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decay amplitudes receive corrections that violate factorization. Phenomenologically the 
most important example of such effects are weak annihilation contributions. Flavor 
non-singlet annihilation effects in B ^ KP decays can be parameterized in terms of 
two parameters /32(MiM2) and /5f(MiM2) 0, Their effects can be incorporated 
by replacing ai ^ ai + P2 and 04 — 0:4 = 0^4 + /3f in the expressions for the decay 
amplitudes. In addition, the term Sp^Sg^u l32{KPs) must be added to the parenthesis 
multiplying fp in (]T3|). 

For the case where the kaon is replaced by a vector meson K*, one makes the fol- 
lowing replacements in the expressions for the amplitudes: fx — > fx*, -^o^^^(O) — >■ 
A^^^*(0), (pxiy) 4>K*{y), and r^* = {2mK*fk*)/{mb fx*), where f^* is the 

(scale-dependent) transverse decay constant of the vector meson. In addition, the rela- 
tion of the parameter 0^4 to the quantities a^^ introduced in [|ll|] changes. We have 



al{PK*) = al{PK*) + r^'al{PK*) in analogy with the relation for al{PK), but 
a^{K*Ps) = a\{K*Ps) — r^^aQ{K*Ps) with a minus sign between the two terms in the 
case of a^{K*Ps). No such sign was present for a^{KPs). Additional changes regarding 
convolutions with twist-3 vector-meson distribution amplitudes, and weak annihilation 



effects, will be detailed in [21 



The predictions obtained using the QCD factorization approach depend on many in- 
put parameters. First there are Standard Model parameters, such as the elements of the 
CKM matrix, quark masses, and the strong coupling constant. Of those, our results are 
by far most sensitive to the strange-quark mass. Specifically, we use \ Vcb\ = 0.041 ±0.002, 
\Vub/Vcb\ = 0.09 ±0.02, 7 = (70 ±20)°, and show results for m, = 100 and 80MeV. Next 
there are hadronic parameters that can, in principle, be determined from experiment, 
such as meson decay constants and transition form factors. In practice, information 
about these quantities often comes from theoretical calculations. The corresponding 
model uncertainties in the form factors have a large impact on our results. We use 
(at g2 = 0) Fff-^'^ = 0.28 ± 0.05, F^^^ = 0.34 ± 0.05, A^^^* = 0.45 ± 0.10, and 
JF^^^'*' = 0.3 ± 0.5 for the "form factors" related to the soft spectator interactions 
discussed in Section Finally, we need predictions for a variety of light-cone distribu- 
tion amplitudes, which we parameterize by the first two Gegenbauer coefficients in their 
moment expansion. Experimental information can at best provide indirect constraints 
on these parameters. Fortunately, it turns out that the sensitivity of our predictions to 
the Gegenbauer coefficients is small, the only exception being the first inverse moment 
of the 5-meson distribution amplitude, for which we take = (0.35 ± 0.15) GeV. (A 
complete list of all input parameters can be found in |TT|, |2T| . 



An important input to our analysis are the form factors Fq~^^^ ' (0) defined in terms 
of the B 7]^'^ matrix elements of the vector current q'y^h. Unfortunately, little is known 
about these quantities. It is not justified to assume that these form factors are related 
in a simple way (modulo SU(3) violations) to the B ^ tt 01 B ^ K transition form 
factors. The reason is that there is a new mechanism available for the case of a flavor- 
singlet meson, where the B meson is annihilated by the current and the singlet meson 
is produced via the emission of two gluons. This possibility was already mentioned in 
Section ^.3|, where it was argued that the two-gluon emission from the light spectator 
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B 



Figure 5: Leading two-gluon contribution to the B 
dot represents the insertion of the current. 



form factor. The 



quark in the B meson gives a leading-power contribution to the B t]^'^ form factor. 
The corresponding diagram, shown in Figure ^ is proportional to the flavor- singlet decay 
constant (v^/p + fp)/V^- We thus adopt the following parameterization for the form 
factors [P = 1] OT rj') : 



Fo«-^(0) = F,| + F,^^i?^±^, (40) 

where Fi and F2 both scale like (A/mfe)^/^ in the heavy-quark limit, and in the FKS 
scheme we expect that Fi ^ F(f^^(0). Given that the ratio iV2 fp + fp)/{V^fn) is 
about 1.2 for P = rj' and 0.2 for P = 77, we expect that the extra contribution might 
be significant for the B ^ rj' form factor, whereas its effect on the B ^ rj form factor 
should be small. In our numerical analysis we set F\ = F(f^'^(0) and take F2 = or 
(somewhat arbitrarily) 0.1. 

Before presenting detailed numerical results for the decay amplitudes it is useful to 
consider the magnitudes of certain flavor topologies in the various decays. Keeping only 
the leading (color-allowed and not CKM-suppressed) electroweak penguin contributions 
for simplicity, we write 

-zA{B- i^-r/W) = {T + C + Pu + S) + A^^) (C, + + Pew + S) , 
-a(P° ^ i^%W) = Ai^) (C + P„ + 5) + A(^) (C, + P, + Pew + S) , 

(41) 

-tA{B- ^ ir-7r°) = Al^) (T + C + P„) + A(^) (P, + Pew) , 



-U(P° ^ K^n') = Ai^) (C - P„) + AW (-Pe + Pi 



EWJ 



The tree amplitude T contains terms proportional to ai, the color-suppressed tree am- 
phtude C contains terms proportional to 0^2, and the penguin amplitudes Pp contain 
terms proportional to and also include the dominant non-singlet penguin annihilation 
contributions. The singlet amplitude S contains terms proportional to , which within 
our approximations does not depend on the flavor label p. The soft spectator-scattering 
contribution as well as flavor-singlet weak annihilation are part of S. The electroweak 
penguin amplitude Pew contains the terms proportional to 0^3 ew small contri- 

bution proportional to a2{KPc) is denoted by Cc- The various topological amplitudes 
depend on the nature of the final-state particles and should thus be labeled T{KP) with 



17 



Table 3: Predictions for the dominant flavor topologies (in units of 10^^ GeV). 
For a given final state the first line refers to F2 = 0, the second to F2 = 0.1. 
Errors are explained in the text. 



Mode 


e^^ X^u^ T 


X['^ P. 


xlf^s 


- , F2 = 

K, ^^ = 
F2 = 0.1 

K-K^ F2 = 


4.0 + O.li 
6.6 + 0.2i 
4.9 + O.li 
5.3 + O.H 
5.8 + 0.2i 


(-65+i?+i?)^(_9+4+i5). 

(-7611^1?) + (-iott;j^)z 

i-29tltt) + 


(-3+i^i) + (4i+^)i 

(-3li^i) + (4+5+2). 

(-0.5llit°i) + (0.6irvlgi)^ 
(-0.5l|^t°i) + (0.6irrlSi)^ 



- ^2 = 

F2 = 0.1 

F2 = 0.1 

i^*7r0 F2 = 


5.4 + O.li 
8.9 + 0.2i 
6.7 + 0.2i 
7.2 + 0.2i 
7.9 + 0.2i 


{-28l,lt,l) + {-lnt\t)^ 
(~29^]^2il9) + ("-'--2^16)^ 


(4^}^;S) + (5^^^g)^ 

(0.6ti^t°) + (0.8l{jt°)z 

(o.6^i^^°) + (o.8^{j^°)^ 




P = rj, f]', vr", etc.; however, in the isospin limit the amplitudes are the same for and 
K'^ in the final state. A decomposition analogous to ( PI| ) holds for the decays B K*P. 

The two most important amplitudes, the penguin amplitude Pc and the tree am- 
phtude T, are given in Table |^ together with the fiavor-singlet amplitude S. For the 
penguin and singlet amplitudes the second error refers to the uncertainty in the cal- 
culation of weak annihilation effects (more precisely, the uncertainty parameterized by 
the quantity Xa), and the first error contains all other parameter uncertainties added 
in quadrature, excluding an overall normalization uncertainty from the CKM factor Vcb- 
In the case of the penguin amplitude the parameter uncertainties are dominated by the 
strange-quark mass. In contrast, the largest uncertainty for the fiavor-singlet amplitude 
comes from the quantities and Xh (defined in [|n[]), which enter the calculation of the 
spectator-scattering effect. In particular, neither of the three effects related to the gluon 
content of singlet mesons, which we identified and calculated in this paper, constitutes 
a dominant uncertainty if the corresponding parameters lie within our estimates. The 
error on the tree amplitude is not given explicitly in the table. It is typically of order 
20% on the real part and of order 100% on the (small) absorptive part, excluding an 
overall normalization uncertainty from the CKM factor V^b- 

The table shows that with the exception of the final state K*r]' the fiavor-singlet 
amplitude S is not an important contributor to the magnitude of the decay amplitudes. 
In a first approximation the pattern of branching fractions is therefore controlled by the 
penguin amplitude Pc, which is indeed quite different for the various final states. We 
first note that the penguin amplitude for the K*7r^ final state is smaller than for Kn^, 
because the K* meson can be produced from a (sg)s+p current only through rescattering. 
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which occurs at next-to-leading order in QCD factorization. ^ The new element for the 
r]*-'-* final states is a contribution to the penguin amplitude in which the kaon (rather 
than the pion or r]^'^) picks up the spectator quark, while the t]^'^ is produced through 
the strangeness content of its wave function. It follows from (|l^) that 

P,{KP) oc Fo^-^(O) ^ aliPK) + F,^-^{0) rp&HKP,) , (42) 

and the interference of the two terms determines the magnitude of the penguin ampli- 
tude. In the case of Krj' the two terms add constructively to yield —20 — 45 = —65 
(referring to the case F2 = and the real parts in Table 0). Since is negative (see 
Table H), the second amplitude changes sign for the Krj final state, resulting in a partial 
cancellation —24 + 37 = 13 of the penguin amplitude. Our calculation therefore con- 
firms the penguin interference mechanism suggested in |12| as the origin of the large Kr)' 



branching fractions. For vector mesons a^iPK*) becomes smaller because the (sg)s+p 
contribution is suppressed, while al{K*Ps) becomes smaller and changes sign because 
the sign of the (sg)s+p contribution is reversed, as discussed at the beginning of this 
section. The two terms in (^2]) therefore give —9 + 20 = 11 for K*ri' and —11 — 17 = —28 
for K*7]. In cases where the penguin amplitude is suppressed by cancellations there is 
a sizeable interfering tree amplitude, which leads to the possibility of large direct CP 
asymmetries. 

Our predictions for the various CP-averaged branching fractions and direct CP asym- 
metries0 are collected in Tables ^ and ^ Since their variations under changes of the CKM 
parameters \Vcb\, \Vub\ and 7 within their respective error ranges are small compared to 
other uncertainties, we have fixed these parameters to the central values specified above 
rather than including them into the error estimates. As in Table ^ the second error 
comes from weak annihilation and the first error from other hadronic parameters. In 
contrast to Table ^, however, we have made the dependence on the strange quark mass 
explicit by showing results for = 100 MeV (default) and = 80 MeV, and do not 
include into the hadronic parameter error. The following observations can be made: 

i) There is a strong enhancement of some branching fractions in QCD factorization 
at next-to-leading order compared to naive factorization (equivalent to QCD fac- 
torization at leading order). For instance, we find BT{K~r]')i^o = compared 
to Br(i^"?7')NLO = 42+}2^ii (for default parameters) and BT{K*'r])i^o = 4l:J com- 
pared to Br(K*-r/)NLO = S.Gll^jt^t" (all in units of 10"^). 



^The smallness of the penguin amplitude for final states containing a vector meson may appear 
problematic in view of the sizeable branching fraction observed for the decay B tt^K*~. The 
complete set of pseudoscalar-vector final states will be discussed in . 

^The sign convention for the CP asymmetry is, contrary to p] , 

Br(g ^ /) - Br(i? ^ /) 
Br(S ^ /) + Br(B ^ /) ' 
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Table 4: Predictions for the CP-averaged branching ratios (in units of 10 
assuming 7 = 70°, \Vcb\ = 0.041 and iKb/K&l = 0.09. The first error is 

due to parameter variations, while the second one shows our estimate of the 
uncertainty due to weak annihilation (see text for explanation). The column 
labeled "default" refers to = 100 MeV and F2 = 0. 



Mode 



Default 



80 MeV 



0.1 



Experiment 



B- 
B- 



■ K-rj 



49+16+27 
-12-11 
+15+26 
-11-11 
1 7+2.0+1.3 
-1.5-0.5 



p-7+21- 
-15- 
9 9+2.7- 
^•^ -2.0- 



1.0 



+1.7+1.1 



q 4+3.2+5.6 
^■^-2.9-2.4 

r A +2.7 +4.5 
^■^ -2.3 -1.9 



c:q+22+41 
"^^-16-17 
-39 
-16 
hl.9 
-0.8 
-1 4+2.4+1.6 
-■-■^-1.7-0.6 
19 +4.3+8.2 
i^-O -3.8-3.5 
o r +3.7+6.8 
-3.1 -2.8 



cfi +19+31 
00-14-13 



56" 



-18+30 
-13-13 
-, 4+1.8+1.1 
-■-■^-1.2-0.5 
r, 7+1.5+0.9 
-0.9-0.4 
q 4+3.2+5.6 
^■^-2.9-2.4 
r q +2.7 +4.5 
^■^ -2.3 -1.9 



72.2 ±5.3 
54.8 ± 10.1 

< 6.9 

< 9.3 
12.7 ± 1.2 
10.2 ± 1.5 



B- 
5° 
B- 
5° 

B- 
B^ 



- K*-r]' 
■ K*-ri 



q r +4.4+4.7 
"^■O -3.7-1.7 



7 7+7.6+8.0 
' ■' -6.7-3.2 



2.5 



+3.8+4.3 
-3.1-1.5 



o R +3.0 +14.0 
°-"J-2.6- 4.4 
Q 7+2.9+14.0 
°-' -2.6- 4.5 
q 9+1.2+4.0 
'-••^-1.1-1.3 
r, 7+0.6+2.4 
•J-' -0.5-0.6 



+6.8+7.4 
-5.8-2.9 
+4.8+1 
-4.2- 



6.3 

-iq 0+4.8+19.8 
l-O.O _4 9 _ g 7 



9 7+3.5+3.9 
-2.6-1.3 
-1 9+2.7+3.2 
-'-■^-1.8-0.9 



13.9 
3.3 



+4.6- 
-4.1 - 



19.5 
6.7 



+1.3+4.8 
'-1.2-1.5 
r, 7+0.6+3.0 
•J-' -0.5-0.6 



+3.1+14.3 
-2.7- 4.6 

Q 9+3.0+14.2 
^•^-2.7- 4.7 
+1.2+4.0 
-1.1-1.3 



9.1 



3.2 



r, 7+0.6+2.4 
-0.5-0.6 



< 35 

< 13 
26.5 ±6.1 
16.4 ±3.0 

< 31 

< 3.6 



) While the calculation reproduces well the qualitative pattern of the known branch- 
ing fractions, we find large uncertainties, principally from weak annihilation and 
the error on the strange quark mass, but also due to the unknown gluon contri- 
bution to the S — > 77^') form factors, parameterized by F2. A smaller value of 
helps to bring the calculation into better agreement with the data. Where cancel- 
lations occur in the penguin amplitude, the uncertainty from weak annihilation is 
amplified. In view of this the significance of the B~ K*"!] branching fraction, 
which comes out too small with default parameters, remains unclear. 

) Given that the K~ri' and K^rj' final states differ only by a CKM-suppressed tree 
amplitude, we are unable to explain significantly different branching fractions for 
these two modes, if 7 = (70 ± 20)°. We therefore suspect that the apparently 
different experimental results should converge. Similar remarks apply to the K*~r] 
and K*°rj final states. 

) Even though direct CP asymmetries are generically predicted to be small in QCD 
factorization (since the strong-interaction phases of the amplitudes vanish in the 
heavy-quark limit), appreciable asymmetries can arise in cases where there is de- 
structive interference between the dominant amplitudes. We indeed find potentially 
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Table 5: Predictions for the direct CP asymmetries (in percent), assuming 
7 = 70° and \ Vub/Vcb\ = 0.09, with error estimates as in Table ^. Experimental 
results are from p], |35| . 



Mode 


Default 


m, = 80 MeV 


F2 = 0.1 


Experiment 


B- 

B'^ 

B- 

B^^ 

B- 

B^^ 


K-T]' 
K-T] 


9+1+4 
^-1-4 
9+1+1 

99+16+24 
-^^-23-20 

in 

7+3+9 
' -3-10 
/I +3+2 
^ -3-2 


9+1+3 
^ -1 -3 

9+1+1 
^ -1 -1 

1 o +15 +23 
-LO _22-19 
o + 8 +6 
^° -14-5 
7+3+9 
' -3-9 
q+3+2 
'-'-3-2 


q+1+4 

'-'-1-4 

9+1+0 

^-1-1 
9(-+21+32 
^"-28-25 

-1 +1 2 +9 

-'-'-'-31-8 
7+3+9 

' -3-10 
A +3+2 
^-3-2 


-4 ±6 

-10 ±8 
3 ±37 


B- ^ 

B'^ 

B- 

B'^ 

B- 

B^^ 


K*-r]' 
K*-r] 


1 o +25 +30 
-'-°-36-25 
7+8+4 
' -15-4 
q+3+25 
'-'-3-24 
A +3+3 
^-3-3 
Q+10+38 
^-10-41 
1 q +15 +24 
-'-'-'-16-28 


-r n+12+16 
-'-'-' -15-14 
_4+4+2 
^-6-2 
/I +2+17 
^-2-17 

q+2+2 

'-'-2-2 

o +9 +42 

_9 _44 

1 9+14+28 
-'-^ -15-31 


q4+48+60 

'-'^-57-41 

1 9+15+ 9 

-'-^-39-12 

q+3+25 

'-'-4-24 

A +3+3 

^-3-3 
Q+10+38 
-"-10-41 
1 q +15 +24 
-'-'-'-16-28 


-5 ±28 
17 ±26 



large direct CP asymmetries for all those modes in Table ^ whose branching frac- 
tion is of order few times 10~® or less. The uncertainties in these predictions are 
very large, typically at least a factor of 2. 

Some of the dominant theoretical uncertainties affecting our predictions cancel in 
ratios of branching fractions. The most interesting ratios are collected in Table ^ and 
compared with the experimental data. For these predictions we perform a complete 
error analysis, quoting (in this order) the uncertainties due to input parameter variations 
(other than m^), weak annihilation, the strange-quark mass, and CKM parameters. We 
find that the ratios are much less sensitive to weak annihilation and the uncertainty in the 
strange-quark mass than the individual branching fractions, facilitating the comparison 
with experiment. In some cases there are, however, still sizeable uncertainties due to 
other input parameter variations, especially those of the heavy-to-light form factors. 

A more detailed verification of the underlying mechanism of penguin interference will 
become possible once the branching fractions for which currently only upper limits exist 
have been measured. We note that we have also computed the corresponding branching 
fractions for AS" = decays, in which the final state contains a pion instead of a kaon 
and a p meson instead of K*. We find that these decays are tree-dominated. As a 
consequence the cancellations among penguin amplitudes are of lesser importance for 
the overall branching fractions, and no drastic signatures as in the case of strangeness- 
changing decays emerge. 
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Table 6: Predictions for some ratios of CP-averaged branching fraction, in- 
cluding error estimates for the variation of all input parameters (see text for 
explanation). 



Ratio 


F2 = 


F2 = 0.1 


Experiment 




A A +2.2+0.2+0.3+0.5 
^•^-1.4-0.2-0.3-0.4 
cr +10.0+0.7+0.1+2.8 
4,7-0.6-0.1-1.7 
n no +0.01 +0.01 +0.00 +0.04 
U-yo -0.02 -0.01 -0.01 -0.04 


a n +2.7+0.2 +0.4+0.5 
'-'•'-'-1.8-0.2-0.3-0.5 
7 n +12.7+1.0+0.4+3.7 
'■'-'- 6.1-1.0-0.3-2.2 

1 nn +0-01 +001 +0-01 +0.05 

^■^^ -0.02 -0.01 -0.01 -0.06 


5.7 ±0.7 
> 1.5 
0.76 ±0.15 


K*-r]'/K*''n^ 
K*^r]/K*-r] 


1 1 +1.3+0.3+1.8+0.7 
^■^ -1.1-0.2-0.6-0.4 
9 7+1.1+0.5+2.0+0.6 
^- ' -0.8-0.5-1.0-0.4 
1 fil +0.04+0.00+0.00+0.17 
-0.05 -0.03 -0.01 -0.14 


n 0+1.0+0.3+1.5+0.6 
'-'•° -0.8 -0.1 -0.4 -0.4 
9 0+1.1+0.4+2.1+0.6 
-^•o -0.8 -0.5 -1.0 -0.4 

1 ri9+0.04+0.00+0.00+0.18 
^■^^ -0.05 -0.03 -0.01 -0.15 


> 0.5 
0.62 ±0.18 



5 Summary 

Motivated by the observation of the large branching fractions for the decays B — > Krj' 
and the distinctive pattern of other decay modes with rj or r]' mesons in the final state, 
we have computed the flavor- singlet decay amplitude in rare hadronic B decays using the 
framework of QCD factorization. We have considered three effects that are specific to 
singlet mesons due to their gluon content, and which occur at leading power in the heavy- 
quark expansion: the b —>■ sgg amplitude, equivalent at leading order to an "intrinsic 
charm" decay constant; spectator scattering involving two gluons; and weak annihilation, 
where the leading-power contribution can be interpreted as a gluon contribution to the 
B r]^'^ form factors, and the remainder can be estimated in terms of the leading-twist 
two-gluon hght-cone distribution amplitude of the r]^'\ A conceptually new result is 
that the spectator scattering is soft at leading order in the heavy-quark expansion, and 
hence appears to break factorization. Nevertheless, the long-distance contributions can 
be parameterized by a non-local B K'^*^ form factor, which is insensitive to the details 
of the singlet-meson wave function. In this sense factorization still holds; however, the 
factorization formula has to be amended for singlet-meson final states by an additional 
term involving this non-local form factor. 

Our numerical results indicate that the fiavor-singlet decay amplitude is not a key 
factor in explaining the pattern of the observed branching fractions, at least if the pa- 
rameters that enter this amplitude (such as the non-local form factor mentioned above) 
do not exceed substantially our estimates. Rather we find (as has been suggested quali- 
tatively in |T^) that the constructive or destructive interference of penguin amplitudes. 



where either the K^*^ meson or the rj^'^ meson picks up the spectator quark, is responsible 
for this pattern. The structure of these cancellations is visible already in the naive factor- 
ization approach [|l^]. The important improvement from QCD factorization comes from 
the possibility to compute radiative corrections to the penguin amplitude, which turn 
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out to be large enough to bring the predicted branching fractions into reasonable agree- 
ment with the data, at least for certain choices of the input parameters. The hadronic 
parameter uncertainties remain however large. We may therefore conclude that while it 
appears unlikely that one can obtain an accurate description of final states with singlet 
mesons from first principles, the results of our analysis clearly support the relevance of 
factorization to this class of charmless hadronic decays. 
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